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We generalize the binary quantum counting algorithm of Lesovik, Suslov, and Blatter [Phys. Rev. 
A 82, 012316 (2010)] to higher counting bases. The algorithm makes use of qubits, qutrits, and 
qudits to count numbers in a base 2, base 3, or base d representation. In operating the algorithm, 
the number n < N = d K is read into a _K"-qudit register through its interaction with a stream of 
n particles passing in a nearby wire; this step corresponds to a quantum Fourier transformation 
from the Hilbert space of particles to the Hilbert space of qudit states. An inverse quantum Fourier 
transformation provides the number n in the base d representation; the inverse transformation is 
fully quantum at the level of individual qudits, while a simpler semi-classical version can be used on 
the level of qudit registers. Combining registers of qubits, qutrits, and qudits, where d is a prime 
number, with a simpler single-shot measurement allows to find the powers of 2, 3, and other primes 
d in the number n. We show, that the counting task naturally leads to the shift operation and 
an algorithm based on the quantum Fourier transformation. We discuss possible implementations 
of the algorithm using quantum spin-d systems, d-well systems, and their emulation with spin- 
1/2 or double- well systems. We establish the analogy between our counting algorithm and the 
phase estimation algorithm and make use of the latter's performance analysis in stabilizing our 
scheme. Applications embrace a quantum metrological scheme to measure a voltage (analog to 
digital converter) and a simple procedure to entangle multi-particle states. 



PACS numbers: 03.67.Ac 03.67.Bg 73.23.-b 

I. INTRODUCTION 

The representation of an integer number n and its de- 
composition into prime factors are basic mathematical 
operations. Quantum mechanics offers a new perspective 
on these tasks, as well known since the seminal work of 
Peter Shoi^ on the efficient factorization of large num- 
bers, with a drastic impact on the security of codes. 
But even the counting of small numbers and their fac- 
torization may prove useful, e.g., in the manipulation of 
physical number states (and superpositions thereof) or 
in the entanglement of flying qubits 2 -"^. A quantum al- 
gorithm to count n < N = 2 K particles propagating in 
a wire using an array of qubits (a iv-qubit register) has 
been proposed recently 4 ; a very similar scheme has been 
proposed by D'Helon and MilburrP in order to find the 
number state distribution of a vibrational excitation in 
a system of trapped laser-cooled ions. Besides providing 
the number n in a binary form (base 2 counting), a sim- 
plified version of this algorithm tests for the divisibility 
of n by 2 k for a given k < K and thus provides the power 
of 2 in the factorization of n. In the present article, we 
generalize this algorithm to perform a base d counting 
and a test for the factor d k in the decomposition of n. 

In order to accomplish this goal, we make use of a mini- 
mal formulation of the counting task in terms of the prob- 
lem of distinguishing between different known quantum 
states in a single-shot measurement. This reduction to 
a few very basic elements naturally connects the count- 
ing task with the quantum Fourier transformation and 
provides us with a constructive scheme for the setup of a 



(non-demolition) quantum counting algorithm. We study 
various possible (hardware) implementations of this algo- 
rithm, paying special attention to the case of a ternary 
(base 3) counting system involving qutrits as elementary 
counting devices. We establish the relation between our 
quantum counting algorithm and the phase estimation 
algorithm and discuss several applications. 

Our counting algorithm is insp ired by the problem 
of counting in mesoscopic systems^' and in quantum 
optics^. A straightforward setup counting particles in 
a non-invasive manner requires of the order of N 2 indi- 
vidual counting elements, see below. A more sophisti- 
cated setup providing a unary counting scheme reduces 
this effort to an order-TV process, and a drastic further 
reduction to a (log d N) 2 scaling can be achieved when 
going over to a counting base d. 

The base- 2 algorithm proposed in Ref. H] involves two 
steps, a first (analog) one where a finite train of n charged 
particles traversing a quantum wire, see Fig. [T] is cou- 
pled to an array of K nearby qubits, thereby rotating 
the states of the qubits in a prescribed manner. In more 
abstract terms, this step corresponds to reading a num- 
ber state |n)* <= in the Hilbert space of number states 
of particles into a state |^ n ) Q of a if-qubit register in the 
Hilbert space H Q of the K qubits. More generally, a su- 
perposition of number states gets entangled with the K- 
qubit register (not copied, as demanded by the no-cloning 
theorem) during the counting process. In a second step, 
the qubit state in H Q is manipulated and read out, pro- 
viding either the maximal power of 2 contained in n or 
the number n in binary form, depending on the readout 
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FIG. 1: Schematic representation of the quantum count- 
ing algorithm (shown is the case of base-two counting with 
qubits). A particle number state |n)* is fed into a quantum 
wire to undergo quantum counting. The interaction between 
the charged particles and the qubits rotates the spin/qubit 
states, thereby generating the Fourier transformation F tak- 
ing the initial state F(|0)q) into the state F(|n)o). A suitable 
manipulation and readout either provides the number state 
|n)cj (this is an inverse quantum Fourier transformation with 
a subsequent single shot measurement in the computational 
basis or a semi-classical quantum Fourier transformation in- 
volving a sequential measurement) or the maximal factor 2 K 
in n (this readout involves 'qubit-rotations' followed by a sin- 
gle shot measurement). 

algorithm. 

The two steps correspond to two quantum Fourier 
transformations of special kind: expressed in the usual 
(cf. Ref. |S) computational basis |0) Q = |00 . . . 0) Q , |1) = 
|00 . . . . . . , \2 K - 1) Q = |11 . . . 1) G H Q of the K- 

qubit register, the passage of the n particles in the quan- 
tum wire transforms the initial state 

i K — l 

l*o) Q = F(|0) Q ) = ^ ]T \k) Q , (1) 

^ fc=0 

i.e., the lowest Fourier harmonic of the -ftT-qubit register, 
into the n-th harmonic 

2 K -1 

l*»>« = F(|n) Q ) = 7^ E exp(2mnk/2 K )\k) Q , (2) 

the quantum Fourier transform of the state \n) Q = 
\n\ ri2 • ■ ■ tik)q, where n is written in binary form, n — 
ni n 2 ...n K = n 1 2 K - 1 +n 2 2 K - 2 + ---+n K 2°, cf. Ref. El 
Depending on the desired information (precise count- 
ing or divisibility) and the available hardware (an oper- 
ating quantum computer or a set of qubits) , the readout 
of the qubit register can be done in various ways. The 
most efficient version for this second step of the algorithm 
is the application of a second (inverse) quantum Fourier 
transformation F _1 on the qubit state F(|n) Q ), taking it 
back into the state \n) Q . A simultaneous measurement of 
the K qubits then provides the particle number n in bi- 
nary form; in case of a superposition of particle number 
states, the measurement will execute a projection onto 



one of them. Alternatively, the number state may be 
used in a further computation. 

However, performing an inverse quantum Fourier 
transformation on the state F(|n) Q ) requires a quantum 
computer (or at least those qubit operations required 
in the execution of a quantum Fourier transformation). 
Instead, we can make use of a procedure which is ba- 
sically identical to the semi-classical Fourier transform 
suggested by Griffiths and MiP^I, a conditional measure- 
ment algorithm involving a sequential readout, where the 
reading of the j-th qubit depends on the results of the 
previous j — 1 measurements. This measurement still pro- 
vides the full information on a pure number state; if the 
initial state is a superposition of number states, the semi- 
classical algorithm will project the state (upon sequential 
measurements of the qubits) to one of its components. 

A specially efficient readout is available if we are in- 
terested in the power of 2 contained in n rather than n 
itself. In this case, a simultaneous (rather then condi- 
tional) readout algorithm can be applied directly to the 
state F(|n) Q ); this algorithm then provides a divisibility 
check of n by 2 , k < K. 

The present article deals with the generalization of this 
algorithm. The most obvious task to generalize is the 
determination of other powers of primes in n, i.e., to 
find the factorization of n, n = 2 k2 3 ks 5 k5 .... This can 
be achieved by going over to generalized qubits, three- 
level systems or qutrits, d-level systems or qudits, etc. 
Equivalently, this corresponds to changing the represen- 
tation of the number n from binary (base 2) to ternary 
(base 3), quinary (base 5), etc. Again, the two-step al- 
gorithm first transfers the information from the physical 
number state (n)^ into the computational X-qudit regis- 
ter (— > F(|n) Q )) through a particle-qudit interaction and 
then extracts the information in the qudit register via an 
inverse Fourier transform (— > \n) Q ). This readout step 
involves a full quantum transform on the level of each 
single qudit, while a semi-classical transform^ suffices 
to extract the information on the level of the X-qudit 
register. 

In order to carry out the above program, it is very 
helpful to have an abstract understanding of the count- 
ing process. Indeed, a minimal abstract formulation of 
quantum counting in an TV-dimensional Hilbert space (al- 
lowing to distinguish or count at most N objects) natu- 
rally leads us to two types of basis states, the computa- 
tional basis (corresponding to the states \n) Q ), in which 
the result of the counting process is measured, and the 
counting basis \ip n ) Q where the actual counting process is 
done — it turns out that just these two basis-sets are nat- 
urally related by the quantum Fourier transformation, 
|VVi)q = F(|n) Q ). Furthermore, the abstract analysis of 
the counting process provides us with a recipe how the 
algorithm can be physically implemented. 

In the following, we compare the efficiency of various 
quantum counting algorithms (Sec. |ll]) and then briefly 
repeat our previous base 2 quantum algorithm with 
qubits, including the sequential and single-shot readout 
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schemes, see section III We then present in Sec. IV our (0/2) sin(n0) has to be much larger than the uncertainty 



basic analysis of the quantum counting process, provid- 
ing a natural link between quantum counting and the 
quantum Fourier transformation as well as a constructive 
scheme helping us to generalize counting to a base d sys- 
tem. We then proceed with the simplest generalization 
to qutrits in section [Vj the problem of counting in a base- 
3 representation or power counting of 3, and discuss its 
further generalization to qudits. To be specific, we con- 
sider an implementation with a three-level system in the 
form of three quantum dots and also discuss various ideas 
for other hardware implementations of the base-3 algo- 



rithm in section VI a spin-1 system, serving rather as a 
Gedanken experiment for illustration, and two practical 
versions of emulating a qutrit with qubits. In Sec. |VIl[ we 
discuss an interesting correspondence between our count- 
ing algorithm and the phase estimation algorith m 1 9 * 11 * 12 ! 
(no such correspondence is yet known for our divisibil- 
ity check) and apply this insight in the proposal for a 
quantum voltage-detector (an analog-digital converter). 
Another application, a scheme to create multi-particle 
entangled states in a Mach-Zehnder interferometer, is 
discussed in Sec. IVIII1 We summarize and conclude in 
Sec. El 



II. EFFICIENCY OF QUANTUM COUNTING 

In a broader context, the efficiency of our quan- 
tum counting algorithm has to be compared with other 
schemes. E.g., the most straightforward non- invasive way 
of counting the number of (charged) particles flowing in a 
wire (directed along x) is to use a spin-counter polarized 
in the xy-plane and rotating the state by an incremen- 
tal angle <j) = k/N (around the z-axis) upon passage 
of a particle^. This spin rotation is achieved throu gh the 
magnetic field pulse generated by the passing charge^^ 
see also Refs. 8 15 for a related experiment in a quan- 
tum optical setup. The precise correspondence of the 
angle with the number of passed particles then requires 
an accurate measurement, to a precision of <f> = tt/N, 
of the spin-counter's final state polarization — this either 
necessitates a large number M > N 2 /ir 2 of repetitions 
of the counting experiment, or a single-shot readout of 
M identical counters all measuring the passage of the 
particles: Measuring the spin along the j/-axis, the (the- 
oretical) probability to find it pointing upwards is given 
by = (rn i )k/M = cos 2 (ncf> / '2) , where (rw)k denotes 
the average of finding of the M spins pointing up in 
a sequence of k — ¥ oo realizations of the entire exper- 
iment. On the other hand, the one-time measurement 
mj„ provides the experimental result P^ = m^/M, from 
which we can find the number n — (2/</>) arccos[(P T J l ) 1 / 2 ]. 
As this procedure is a statistical one, we have to deter- 
mine how many spins (measurements) M are needed to 
predict the particle number n with certainty. We then 
require that the difference in probability (we assume 
N > n > 1) = |pt( n + 1) - pt( n )| ~ \Q n pt\ = 



[((W) 2 ),] 1 / 2 



[((mt 



)fe) 2 )fc] 1 ^ 2 in the measure- 



ment, SP^ > [((<5m 1 ') 2 ) fc ] l / 2 /Af. Given the binomial 
statistics of the measurement process [the values f an d 
I are measured with probabilities P^ and (1 — P^)], we 
obtain ((8m^) 2 )k = P^(l — P') M and combining these 
results, we find that 



M > l/4> 2 > N 2 /n 2 > 1 



(3) 



spins are needed in order to accurately measure the par- 
ticle number n < N. 



III. QUBITS: COUNTING POWERS OF 2 

We first provide a more detailed discussion of the base- 
2 counting algorithm with the setup in Fig. [TJ where the 
n < N = 2 K particles to be counted flow in a quantum 
wire along x. Single elect ron pulses can be generated by 
appropriate voltage pulse d 6 * 16 * or through injection from a 
quantum dol^and the counters are conveniently thought 
of as individual spins, cf. Refs. l6|13)14l We use spin states 
polarized along the z-axis as our computational basis, 
|t) |0) and \i) <H- 

Preparation: Initially, the K spins or qubits (we use 
these terms synonymously, cf. Ref. 18) are polarized 
along the positive y-axis, i.e., the initial states read 
I + V)j = [| t), + i\ l);]A/2, j = l,...,K. Iden- 
tifying | |)j <->• |0)j and | l)j O — the prod- 
uct state |* )q = ]ljLi[(|0)i + |l)i)/V^] is identical 
with the equally weighted sum of K-qubit register states 



|o>„,|i; 



■1 
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I* 



0/Q 



= (l/V2*)£ fc |fc) Q . This 



state then coincides with the lowest harmonic in the 
Fourier transformed computational basis: indeed, the 
quantum Fourier transform takes a state 



(4) 



into the state 



F(\X) Q ) = J2yk\k) Q = \Y) Q (5) 
k 



with 



y k = {l/VN)J2xi exp(2wilk/N). 



(6) 



The initial state \^o)q of the counting setup then is given 
by xi = Sio and hence |* )q = f (I°)q)- 

Counting and Fourier transformation: Assuming a 
transverse coupling between the charged particle and 
the spin, the passage of a particle rotates the spins in 
the x-y plane. The couplings of the spin counters to 
the wire are chosen such that the j-th spin is rotated 
(anti-clockwise) by the amount <j)j = 2n/2 3 (a rotation 
by U z ((f)j) = exp(—i(j>jcr z /2) with a z a Pauli matrix). 
The passage of n particles then rotates the j-th spin 
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by the amount n<pj and thus it ends up in the state 
[|t)j + iexp(2mn/^)\X)j]/V2, where we have dropped 
the overall phase exp(— ni n/2- 7 ). Again, we identify 
| f)j «-> |0)j and | l)j — *| l)j and use the binary 
representation of integers k — k\l K ~ x + • • • + fc_f<-2 = 
k\ k-2 . . .fcif to rewrite the product |^ n )Q = + 
exp(27rin/2 3 ')|l) J -/2] over qubit states as a sum over reg- 
ister states, 



7 6 5 4 3 2 1 



K 



l*n) Q = II ,V "-' 



|0)j +cxp(27rin2^-V 2K )l 1 )j 



V2 



(7) 



V" e 2^in(J2^ 1 k l 2 K - , )/2 K 
^fa,.., 

1 53e^- fc ^|fc>g 



1*1 



fci,...,fejf =0,1 
2 r 



From the comparison with the Fourier transform Eqs. 
{5} and @ we find that xi = <5/„, and hence the pas- 
sage of the n particles (the counting process) takes the 
state F(|0) Q ) into the n-th harmonic F(|n) Q ). We call the 
states \n) Q the computational basis and the transforms 
\^n)q = F(|n) Q ) define the counting basis. The fact that 
the qubits in the if -qubit register resid e in a product 
state and hence remain unent angle d 1 ° * 1 2 1 is a crucial el- 
ement of our algorithm and in fact decisive for the next 
step, the readout of the result with the help of a semi- 
classical quantum Fourier (back) transformation. 

Readout and inverse Fourier transformation: Let us 
then turn to the second step, the semi-classical quantum 
Fourier (back) transformation which provides us with the 
desired result, the binary representation of n. We start 
with the measurement of the first spin: since this has 
been rotated by the angle nn, we measure it along the 
y-axis. If we find it pointing upward, the number's parity 
is even and we store a '0' in the first position uk of the 
binary number; in case we find it pointing downward, the 
parity is odd and we store the digit '1'. Besides providing 
the number's first binary digit, the parity, the outcome 
of the measurement also tells us whether the second spin 
(rotated by nn/2) is directed along the y-axis (even par- 
ity) or along the x-axis. This information allows us to 
measure the second spin along the correct axis; for an 
even n, we measure the spin along the y-axis (and store a 
(1) in riK—i if the spin is pointing up (down)), while for 
an odd-parity n, we measure the spin along the — re-axis. 
The iteration of the readout algorithm is straightforward: 
the j-th spin is measured along the direction m,j—i<f>j with 
the integer rrij-\ = n K-j+2 ■ ■ ■ nK-inx corresponding to 
the binary number encoded in the j — 1 previous mea- 
surements. The j-th position in the binary register then 
assumes a value n,K-j+i = or nx-j+i = 1, depending 
on the measurement result, for a spin pointing along 
the axis and 1 for a spin pointing opposite. This sequen- 
tial measurement algorithm provides us with the binary 
representation of n. 

In order to render the algorithm more efficient, rather 
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FIG. 2: Counting of n — 7 particles by K = 3 qubits. Af- 
ter passage of the particles, the initial state with all spins 
pointing along the y-axis and corresponding to the lowest 
Fourier harmonic is transformed to the final state, the n-th 
Fourier harmonic with properly rotated spins. The indivdi- 
ual counter qubits are then measured (M) and the binary 
digits nj, j = 3,2,1 are determined in reverse order. After 
each measurement, the remaining qubits are rotated (R) to 
undo the counting of the corresponding moduli: Following 
the measurement of the first qubit, the rotations by the last 
(odd numbered) particle are undone. After measurement of 
the second qubit the rotation by the last odd-numbered pair 
is undone. 



than rotating the axis of measurement, the spins are 
rotated backwards by the corresponding angles. These 
rotations by —rrij—ify are conveniently done incremen- 
tally: after measurement of the j-th spin with outcome 
'0' or '1', all spins J > j are rotated by —riK-j+i^^^J- 
These rotations undo the action of odd-numbered groups 
of particles: The first rotation by —2wriK 1/2 J acting on 
the qubits J > 1 compensates for the last odd-numbered 
particle. The second rotation by — 27rnx-i2/2 J act- 
ing on the qubits J > 2 compensates for the last odd- 
numbered pair. The third rotation by — 27rn^_24/2 l/ 
acting on the qubits J > 3 compensates for the last odd- 
numbered quartet, etc. These rotations make the next 
spin to be measured point up or down, since the action 
of those particles giving an intermediate result has been 
subtracted. The entire algorithm is illustrated in Fig. [2] 
for the case of n = 7 particles counted by K = 3 qubits. 
Formally, the availability of a sequential readout algo- 
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rithm can be derived from a suitable representation of 
the product state Eq. (Frk: the fraction n/2 J in the phase 
exp(27ri n/2 J ) of the j-th qubit has to be known only 
modulo 1 and making use of the relation 



n 



Dd(l) 



= 0.n K - 



3+1 



,n K 



(8) 



riK-i+i2 j 1 



2J 



the j-th qubit state can be written in the form [\0}j + 
exp(27ri O.riK-j+i ■ ■ ■ n x)\l) j]/V%, where we make use of 
the binary representation of fractions 0.nin 2 ■ ■ ■ Uk = 

m/2 + n 2 /4 H h n K /2 K , see Ref. M The final state 

after passage of the particles can then be written in the 
form 



K 

i*n>« = n 



|0)j + cxp(27ri 0.rtK-, + i ■ ■ .n K )\l) 3 
V2 



(9) 



This representation demonstrates that the state of the 
first qubit j = 1 involves only the smallest digit uk of 
the seeked number n, the second one involves the fraction 
O.nx-inji, and so on. Hence each qubit state requires 
knowledge of the states of previous qubits and its mea- 
surement adds one digit more to the binary representa- 
tion of n. 

Divisibility by 2 k : A variant of the above counting al- 
gorithm provides a test for the divisibility of n by powers 
of two. Consider the state of the first j = 1, 2, . . . , k < K 
spins after the passage of n particles. If the number n 
contains the factor 2 fc , then the k qubits will all point 
along the positive y-axis (and, for k < K, the k + 1-th 
qubit will point down; alternatively, if n = 2 K then all 
spins in the register have returned to their initial state 
pointing along y). A single-shot measurement of the K- 
qubit register along the y-axis thus provides the (maxi- 
mal) factor 2 k in n. The formal proof of this statement 
is given in Ref. HI 



IV. DISTINGUISHABILITY AND QUANTUM 
COUNTING 

We now reduce the problem of quantum counting to 
the task of distinguishing between quantum states. This 
reduction will quite naturally lead us to the definition 
of two basis sets, one serving the counting process itself 
by admitting a simple manipulation of phases during the 
counting step (the computational basis) and the other 
keeping track of the counting (counting basis); the two 
are related by the operation of quantum Fourier trans- 
formation. 

We start from the premise that quantum counting cor- 
responds to the process of associating distinct states of an 
auxiliary quantum system (the counter) to the size (car- 
dinality) of a physical state. We assume that we want to 
count at most N objects, hence our auxiliary quantum 



system shall count the objects modulo TV. During the 
counting process, the initial state |\l/o) of the counter is, 
upon passage of n objects, transformed to the final state 
\^ n ) (we can safely drop the index Q in this section). We 
define the unitary operation Ci to describe the passage 
of one particle, 



= Ci|* > 



and a simple iteration produces the state 



|¥„> = CJ|* ) = Q|*c 



(10) 



(11) 



upon passage of n particles. We now require that we 
can distinguish between the states |\& n ) and |^ ) m a 
single-shot measurement, implying that the states should 
be orthogonal, (^nl^o) = 0. So far, we only require 
that we can distinguish between 'no particles' associated 
with the state an d a state with 'some particles 
with < n < N and associated with the state \^ n ), 
without being able to decide between different number 
states with different n's. It turns out that a setup solving 
this reduced task is also able to distinguish between the 
different particle number states \^ n )- Indeed, using Eq. 
(fill and the fact that (*„|* > = for all < n < N, 



we find that (we choose < I < n < N) 



(*o|C|C„|* ) = 

(*o|*n-J> = 



(*o|C„-i|* > (12) 



and hence the states \^> n ) are all orthogonal and distin- 
guishable. Finally, our wish to count modulo N requires 
cyclicity, i.e., 



C N \^ ) = expi9|* ) 



(13) 



if the dimension of our auxiliary counter system is given 
by N, then the cyclicity Eq. (13) follows automatically: 
Applying Ci iteratively to the counting states l^n), one 
finds that all states |^ n _|_i) are orthogonal to the previ- 
ous states |* ; ), (*j|* n+1 ), 0<l<n<N-l. Once we 
arrive at n = N— 1, the last state completing the Hilbert 
space, the further application of Ci produces a state 
\^n) — CiI^jv-i) which has to be a superposition of the 
previous states, \V N ) = Yln^oi^nl^N^n)- However, 
since all matrix elements (\& n |^jv) = (^nlCi^iv-i) = 
(^n-il^iV-i) vanish for < n < N — 1, we must have 
\^n) c« l^o) and since Ci is unitary, we arrive at the 



result Eq. (13). Otherwise, for a larger dimensionality 



of the auxiliary system, the condition Eq. ( 13 ) has to be 
imposed as a separate requirement. 

In the end, the (minimal) auxiliary counter system is 
described by an A^-dimensional Hilbert space % with or- 
thonormal counting basis |^ n ) £ H, n — 0, . . . , N — 1 
and (if>i\i& n ) = Si n , and a unitary (shift or count- 
ing) operator Ci taking one counting state to the next, 
Cil^n) = l^n+i), and the property of cyclicity, = 
Cjy = exp(iO). In a specific physical implementation, 
the phase O is determined by the dynamical evolution of 
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the system during counting. The states \^ n } keep track 
of the numbers in the counting process, i.e., due to their 
orthogonality they uniquely identify the cardinality of 
the counted set. 

In order to further characterize the properties of our 
auxiliary counting system, we determine the eigenvalues 
and eigenvectors of the counting operator Cx. Expressed 
in the basis {{^ n)}n=o > the latter assumes the form 
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Its eigenvalues and eigenvectors are easily found: the de- 
terminant of Ci — Al is given by (—A) + (— l) N ~ 1 e 
and hence the eigenvalues of Ci are the N roots of 
1 on the unit circle in the complex plane shifted by 
Q/N, A fc = exp{2irik/N + iQ/N), k = 0, 1, . . . , N - 1. 
The associated eigenvector |fc) is given by (^> n \k) — 
exp(-27ri kn/N - i n&/N)/y/N. Note that a phase 6 = 
2irl simply renumbers the eigenvectors and eigenvalues 
by I. 

The eigenstates \k) of the counting operator C x show 
a particularly simple behavior in the counting process — 
they merely pick up a phase, and these phases are dis- 
tributed homogeneously over the unit circle. Hence, ex- 
pressing the (unknown) counting states |^„) through the 
eigenstates \k) of the counting operator Ci, we obtain 

i 

|* n ) = _^J2e 2 " kn / N+lne / N \k) (15) 

= e ln0/N F(\n}). 

Making use of the eigenstates the counting step can 
be implemented in an extremely simple and minimally 
invasive manner: upon passage of a particle, each state 
|fc) shall pick up the additional phase exp[i(27r k + 0) /N] 
and \ty n ) goes over to which is just the action of 

Ci or of counting. Hence, if we choose as our computa- 
tional basis the set of orthonormal eigenstates {\k)}^~Q 
of the shift operator Ci, then the counting process can be 
implemented in a 'soft' way, adding only phases to the 
computational states (note that these are the states in 
which our final projective measurement will be done). 
At the same time, the counting basis, which is inter- 
connected by the shift operator Ci, is made from the 
states \& n ) = exp(i n<d / N)F (\n)) , which is nothing but 
the quantum Fourier transform (up to a phase) of the 
eigenstates |n). In the further general discussion below 
we set = 0; the required transformation eliminating a 
finite will be discussed later. 

At this point one may ask if there arc other count- 
ing algorithms which do not exploit the quantum Fourier 
transformation — the answer is yes, but such alternative 
schemes do not provide the 'soft' counting involving only 
the addition of phases. As an example, consider the setup 




FIG. 3: Sketch for a real space counter in the form of a 3- 
level qutrit operating in the amplitude mode; the interaction 
between the counter particle and the particle in the quantum 
wire is assumed to be attractive (a slightly modified version 
can be found for a repulsive interaction) . Subsequent lowering 
of the barriers separating the quasi-classical states |0), |1), and 
1 2} upon passage of a particle moves the counter particle to 
the left, cf. (a) — > (b) — > (c). The third particle passing drags 
the counter particle all the way back to the initial state, cf. 
(d). 



shown in Fig. |3j a multi-well system (one counter particle 
in a d-well potential landscape) operating in the ampli- 
tude mode. Assuming an attractive interaction^, the 
passing particle lowers or removes the barrier between 
adjacent semi-classical states, allowing the counter par- 
ticle to move between two wells. Here, the semi-classical 
states localized in the individual wells play the role of 
the counting states |^ n )' initializing the counter in the 
(right-most) state \^o), cf- Fig- [§) the first d — 1 parti- 
cles passing will push the counter particle to the left in 
steps of one, |*$ n ) — > |^ n +i), while the d-th particle will 
drag the counter particle back to the right until it ends 
up in the initial state |^o)- The difference of this device 
with the 'soft' counting device discussed above is in the 
choice of the computational basis: rather than selecting 
the Fourier transformed states of \^ n ), which are eigen- 
states of the shift operator Ci and only pick up phases 
during the counting process, here we choose as a compu- 
tational basis the counting basis itself, hence, \n) = \$f n )- 
As a result, rather than adding phases during counting, 
we shift the counter particle in real space. Obviouly, this 
setup is difficult to realize as quite some fine tuning is re- 
quired to generate a clean shift operation; furthermore, 
the shift operation in real space will generate an appre- 
ciable back action on the passing particles. On the other 
hand, such a counter is not supposed to evolve coher- 
ently between counting steps, hence the requirements on 
the coherence time are reduced. 

In abstract terms, the counting process can be illus- 
trated through the counting states \^ n ) arranged in a 
circle with the shift operator Ci transforming one state 
to the next. The goal then is to have for the Fourier 
transformed basis states | k) a set of (measureable) states 
which merely pick up phases when interacting with the 
passing particles; these states then shall form our com- 
putational states. In this basis, the counting process (the 
shift operator Ci) transforms one Fourier mode into the 
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next, cf. Fig. [4^a). 



(a) l¥,) [%) 




FIG. 4: (a) Counting device for N — 8 (corresponding to 
a qudit with d = N = 8). The (semi-classical) states \k), 
k = 0, . . . , N — 1, form the computational basis; their Fourier 
transforms \^ n ) define the counting basis. In the elementary 
counting step, the shift operator Ci transforms one Fourier 
mode |^ n ) into the next \ty n +i). (b) Emulation of the N = 8 
qudit through three qubits, implementing quantum counting 
with d = 2 and K = 3 (left) and corresponding classical 
Abacus (right). 

So far, the abstract computational states \n) have been 
chosen in a trivial way, without any additional structure; 
correspondingly, our counting process is still a simple 
(and inefficient) one, requiring as many different states 
as we have objects to count. The great reduction in 
the (hardware) complexity of the counting process ap- 
pears with the introduction of a counting base. Intro- 
ducing a counting base d (d — 2 (3) for binary (ternary) 
counting), we have N = d K and can reduce the hard- 
ware requirement from N to log d N = K; correspond- 
ingly, the efficiency of the algorithm increases dramat- 
ically from a linear in N complexity to log d N. To il- 
lustrate the case, we consider a simple example N = 8, 
d = 2, and K — 3, counting up to 8 with the help of 
3 qubits. The primitive version involves a loop with 8 
states |0), |1), . . . , 1 7). This can be reduced to three loops 
with two states each, |0) 3 , j = 1, 2, 3, cf. Fig. [4] The 

preparation of the large loop generates the state I'&q ) = 
Sfc=o |k)/\/8, the properly prepared reduced system is 

in the state l^ 3 ') = nf=i[l°)j + \ Thc count- 

ing step on the large loop adds the phase exp(27ri/c/8), 
\k) -> exp(27rifc/8)|fc), to the state \k) in |*[, 8) ). On the 
other hand, the counting step on the small loops adds 



sequentially smaller phases: \u)\ — > ex$(2itiv /2)\v)\, 
\v} 2 exp (27^/4)1^)2, \v) 3 -> exp(27rw/8)|z/) 3 , v = 
0, 1. One easily checks that the state 

l^ 3 ), = [|0)l + |l)l]<8 [|0)2 + |l)2]<g>[|0) 3 + |l)3] 

10 ; V8 

= [|000) + |001) + |010) + |011) (16) 

+ |100) + |101) + |110) + |111>]/V8 

(the three entries refer to the qubits 1, 2, and 3, e.g., 
|011) = |0il 2 l 3 ) = |0> 1 1 1) 2 1 1>3) transforms to the state 

|^f 3) ) = [|000> + e" /4 |001) + e I7r/2 |010) (17) 
+ e 3i7r / 4 |011) + e 47r |100) + e 5 " /4 |101) 
+ e 3 " /2 |110) +e 7l " /4 |lll)]/\/8 

upon passage of one particle; this is identical to ) 
upon identifying the state \k) with its binary equivalent 
\k1k2k3). Hence, the complexity of the hardware used 
in the counting task is easily reduced by going over to 
a counting base; in our case, rather than implementing 
an 8-level system, the same counting task can be accom- 
plished with the help of 3 qubits. This step in reduction 
of complexity is nothing but going over to a (quantum) 
Abacus, cf. Fig. [4jb) . Its physical implementation in- 
volves qudits (playing the role of the rows in the classical 
Abacus); performing one cycle in the j-th qudit shifts the 
state of the next qudit by one unit. 

The insight provided by the 'soft counting procedure' 
and the reduction in hardware provided by choosing a 
counting base (base 2, 3, . . . , d) described above gives 
us a recipe how to construct a physical implementa- 
tion of the counting process: Define a quantum me- 
chanical system with = d orthogonal states |/c), 
k = 0, ... ,N — 1, with an identical (trivial) time evo- 
lution; these states form our computational basis and 
often appear in the form of semi-classical (measureable) 
states. The non-trivial time evolution of these states orig- 
inates from their interaction with the particles during 
their passage. The coupling of the states to the par- 
ticles has to be arranged in a way such that the state 
\k) picks up a phase ex.p(2irik / N) upon passage of one 
particle. Defining the interaction Hamiltonian Hi nt , this 
implies an equidistant distribution of the matrix elements 
jl k dt(k\\-\i n t\k); assuming that all levels interact with 
the particles during equal time intervals tk = t c , we can 
conclude that the spectrum of H; nt is equidistant. Fur- 
thermore, returning back to the cyclic phase O in the 
counting operator Ci, we can compensate for a finite 
value via an energy shift in the interaction Hamiltonian 
H int : adding a constant energy Qh/t c to H int , we rede- 
fine the counting operator operator Ci = e~ lB ^ N Ci and 
the counting basis |§„) = e^ i0 "/ JV |*„), < n < N, 
and find that the new condition for the cyclicity reads 
|*jv) = Cil^-i) = e- ie |^) = e- 4e e ie |*o) = |*o> ; 

As a specific example, we can examine the situation 
for the base 2 counting in Sec. |III| here, the Hamilto- 
nian (more precisely, the logarithm of the shift operator, 







FIG. 5: Equidistant spectrum for the Hamilonian Hi nt de- 
scribing K = 3 qubits. 



(l/i)lnCi) describing one count can be written in the 
form 



If 1 " K 7T 



(18) 



where the Pauli matrix o~ ( ^ operates on the j-th qubit, 
resulting in the expected equidistant spectrum, cf. Fig. 
[5j In this example, the energy zero is located in between 
2 eigenstates and hence = tt. Shifting the zero-energy 
point to one of the eigenstates, we then can get rid of 
the cyclic phase 0; e.g., shifting the energy zero to the 
top-most level we identify the state | ftt) with the com- 
putational state |0). 



V. QUTRITS: COUNTING POWERS OF 3 

As a first step towards the generalization to qutrits, 
we reformulate the base 2 counting algorithm in terms 
of manipulations of a particle in a double well potential 
with semi-classical (computational) states |0) and |1), cf. 
Fig. [61(a) , as directly realizable with a double-dot charge 
qubipffiS. 

We assume the double-dot charge qubit to be aligned 
perpendicular to the wire, such that the two wells couple 
differently to the charge of the passing electrons, cf. Fig. 
[6]ja). We work with the quasi-classical states |0) <-> | f) 
and |1) <-» ||) (our computational basis) and consider 
a 'phase mode operation' of the counter, with a large 
barrier separating the quasi-classical states, resulting in 
an exponentially small tunneling amplitude oc A, with 
2A the gap between the two true eigenstates. The qubit 
is manipulated by applying voltage pulses, either to lower 
the barrier between the two wells in order to change the 
amplitude (rotation around the x-axis; e.g., opening a 
gap 2A between the qubit eigenstates during the time 
t = tm/2A moves the state |0) to the state |1)) or to 
disbalance the two wells in order to change the phase 



FIG. 6: Double- and triple-well potentials defining qubits 
(a) and qutrits (b). The qubit is manipulated via voltage 
pulses 14 lowering the barrier between the wells and chang- 
ing the amplitude, or voltage pulses V v disbalancing the wells 
and changing the phase of the qubit's state (and similar for 
qutrits) . 



(rotation around the z-axis; e.g., lifting the right well 
by 5 during the time t = tup/ 8 adds a relative phase 
exp(— iip) to the state |1)). 

The algorithm involves three steps: To prepare the 
qubit, we start from the semi-classical state |0) and 
apply the unitary (amplitude-shift) operator Aio = 
exp(— iHiot/h) with the tunneling Hamiltonian Hio = 
-A(|0)(1| + |1)(0|) and the time t = hn/4A to produce 
the balanced superposition |\&o) = ((|0) +i|l))/\/2, 



Aio = exp(— iH 10 i/7i)|t=fcir/4A 
1 / 1 i 

~~ 7i \ * 1 



(19) 



in a spin language this corresponds to a rotation around 
the x-axis by —n/2 to produce the initial state | + y). 

Next, we let the particles pass the (first) counter qubit. 
Upon passage of one particle, the semi-classical qubit 
states \v) pick up a phase exp(27rw/2). The unitary op- 
erator Ci = exp[wr(0|0)(0| + 1|1)(1|)] then takes |* ) 
into = Ci|*o) = (|0> +ie i7r |l))/v%; the passage of a 
second particle brings the state \^\) back to |\&o), hence 
Ci is cyclic, Cf = 1. This operation corresponds to the 
rotation of the spin by n around the z-axis. The states 
l^o) an d l^i) define the counting basis. 

Finally, we translate the measuring process. The ap- 
plication of the operator A^g 1 = A} transforms l^o) back 
to |0) and to |1) (up to a phase); a simple check that 
the qubit is in the state |0) tells, that the number n of 
passed particles is even. 

The subsequent qubits j > 1 are prepared in the same 
way. Since it is the task of the qubits j > 1 to detect the 
passage of groups of particles, these counters are more 
weakly coupled to the wire. In particular, the operator 
Ci for the second qubit measures the passage of pairs, 
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hence the phase added to the state |1) is ir/2, and sim- 
ilar for the following qubits. Regarding the final step of 
measurement, there are two variants: in order to extract 
the maximal factor of 2 k in n, it is sufficient to apply 
the operation A| q and then check, whether the counter 
resides in the state |0); the first qubit residing in the 
state |1) determines the power k. On the other hand, if 
the goal is to find the cardinality n, then before apply- 
ing A|q the qubits have to be properly prepared through 
a 'rotation around the z-axis' (rather then 'rotating the 
axis of measurement'). These 'rotations' are conveniently 
done incrementally and involve a phase shift of the quasi- 
classical state |1): after measurement of the j-th qubit 
with outcome '0' or '1', the measured value is stored as 
the digit nx-j+i (note the reversal in the sequence) and 
all qubits J > j are given a phase shift — nA'-j+i2-' _1 0J 
on the state After application of A{ , the qubit j + 1 
is measured. 

An interesting subtlety concerns the possible entan- 
glement of the qubits with the passing particles dur- 
ing counting. The spin counter discussed in section III 
is entirely unproblematic regarding this aspect, as the 
gauge interaction leaves the counted particle essentially 
unchanged. On the other hand, a charge qubit acting as a 
counter can become strongly entangled with the counted 
particle: With the qubit charge in a quantum superpo- 
sition, with 1 /2 probability to be closeby the wire (state 
| t)) and 1/2 probability to be further away | \), the par- 
ticle passing by is first decelerated and then reaccelerated 
by the qubit charge. This process happens with different 
strengths depending on the qubit state. As a result, the 
particle wave function may split after passing the qubit, 
with one part (the fast one, |/)) moving ahead of the 
other (the slow one, \s)). The calculation of the prob- 
ability V y to find the counter state pointing along the 
y-dircction is given by the partial trace over the particle 
space; if the two states |/) and |s) are distinguishable, 
(/I s ) = 0, then V y = 1/2 independent of the phase (j> 
picked up by the counter qubit. Hence, it is crucial that 
the counter does not generate a wave function splitting 
when the particle passes by, i.e., the counter only works 
properly if |/) ~ \s) and the final state is essentially 
non- entangled. Using the charge qubit as a measuring 
device, the requirement of weak splitting boils down to 
the condition^ 4> ^ £^f, where </> is the angle of the 
qubit rotation quantifying the qubit-particle interaction, 
k F is the Fermi wave vector in the quantum wire, and £ 
is the width of the wave packet. 

In a more quantitative analysis, we can consider the 
evolution of a (Lorentzian, cf. Ref. [TB|) wave packet 
^(x) with Fourier amplitudes f(k) — y/^ir^ exp[— (k — 
fc F )£]0(fc — fc F ) subject to the scalar field of a charge 
qubit. For simplicity, we assume that the state |0) acts 
with a potential V(x) on the wave function, while the 
state |1), which is further away from the wire, has no 
influence on the particle. Furthermore, we choose the 
potential V(x) such that the qubit state is rotated by 
(/>, i.e., after the passage of the particle, the initial qubit 



state |* ) = (|0) + i\l))/V2 shall be rotated into the 
state |*0) = (e l ^|0) +i\l))/^/2. Assuming a smooth po- 
tential V(x), we determine the asymptotic particle-qubit 
state within a semi-classical approximation and project 
the wave function onto the qubit state l^^,}] tracing over 
the particle degree of freedom, we find the probability 
— 1 — 2 /8/Cp£ 2 to measure the qubit in the state \%!$) 
and the second term provides the probability to observe 
a wrong result. Note that the error probability involves 
the square of the small parameter </>/£fc F . 

The base 3 counting with qutrits in the form of triple 
wells follows the same scheme as the base 2 counting with 
double-dot qubits described above. We start out with one 
qutrit counter, i.e., K = 1, which initially resides in the 
state |0), cf. Fig. [6^b). First, the qutrit is prepared in 
a balanced state with equal weights in each of the three 
semi-classical ground states (the phases may be chosen 
arbirarily). The (amplitude-shift) operators A 10 and A 21 
performing this task generalize the operator A10 above. 
The angle xio = Aiot/H in the operator A10 (we write 
the matrices in the semi-classical basis), 

A 10 = exp[ixio(|0>(l| + |l)<0|)] (20) 

cosxio isinxio 
isinxio cosxio ] , (21) 
1 

is chosen such that 2/3 of the wave function is shifted to 
the state |1), hence tanxio = V2. The matrix A10 then 
assumes the form 



Aio 




1 iV2 
iV2 1 
\/3 



(22) 



The operator A 2 i transfers weight between states |1) and 
|2), 

A21 = exp[»X2i(|l>(2| + |2)<l|)] (23) 
10 
= I COSX21 «sinx2i ] • (24) 
«sinx2i COSX21 

Its task is to take the state Ai O |0) = (|0) + i\/2|l))/\/3 
to the balanced state |* ) = (|0) + e ilfil |1) + e ilfi2 |2))/v / 3 
(with approriate phases ipi,2)- Choosing the time such 
that X21 = A 2 i i//i — it/ 4 we obtain the operator 



A21 




(25) 



and the balanced state l^o) = A21A10IO) reads 

|*o> = A 21 A 10 |0) = (|0) + - |2))/V3. (26) 

We call the combination U p = A 21 A 10 the preparation 
operator. The particular construction of the state |\&o) 
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does not produce the lowest Fourier harmonic \^fh) = 
[|0) + |1) + |2)]/v3; if we insist to work with the canoni- 
cal expression for the Fourier transformation, we have to 
properly redefine the phases of the computational basis 
states, |0) -» |0), |1) ->■ and |2) -> -|2). Other- 

wise, we can actually start our counting basis with any 
Fourier harmonic of the computational basis or even with 
any other balanced state as it is produced in a convenient 
physical preparation step (as we have done above). Also 
note that our preparation operator U p is not acting as a 
Fourier transformation on the other basis states but is a 
simpler operator that only transforms the state |0) into 
a properly balanced state. 

Next, we find the wave functions |^i) and l^a) after 
passage of 1 and 2 particles: Upon passage of one par- 
ticle, the semi-classical qutrit states \v) pick up a phase 
exp(2mv/?>). The unitary operator 

d = exp[(2™/3)(0 |0)(0| + 1 11)(1| + 2 |2)(2|)] (27) 

generates the additional counting states 



|*i> = Ci|tf > 

|*2> = Ci|*l) 



+ i e 2ir</3| 1 ) _ e 4W3| 2 ) 

+ ie 47ri / 3 |l) - e 87ri /3| 2 ) 



(28) 



(29) 



The further application of Ci brings us back to |\&o)j i-e., 
Ci is cyclic, Cf = 1 (i.e., = 0). It is easy to check that 
the set 1^0)1 1^1)1 1^2) forms a new orthonormalized 
basis in the qutrit 's Hilbert space spanned by the semi- 
classical states |0), |1), 1 2). 

We come to the readout step. We define the inverse 
preparation operator U" 1 = [A21A10]™ 1 = AJqA^; using 
the explicit forms in Eqs. (22) and (25) we obtain the 
expression 



a/2 -iV2 -V2' 
IT [ = -4= I -2i 1 -i 
-iy/3 V3 



(30) 



The application of U p 1 takes the state |^o) back to 
|0), while U~ 1 |\E , i i 2) are still superpositions of the semi- 
classical states |1) and |2), 

u; 1 !^) = |0), (31) 
U; 1 !^) = -i(\l)-\2))/V2, (32) 
U- 1 !^) = -i(|l) + |2»/V2. (33) 



The divisibility test by 3 then involves the application of 
Up 1 to the state C^l^o) obtained after passage of the n 
particles and a measurement of the qutrit state; if the 
qutrit resides in the state |0), then n is divisible by 3. 

Finding the modulus of n to the base 3 is slightly more 
involved. We have to rotate the states UZ" 1 1 ^ 1,2) such 
as to recover the original states |1) and |2) of the com- 
putational basis. Using the spin language in the two- 
dimensional space spanned by |1) and |2), this opera- 
tion is achieved by a 'rotation' by ir/2 around the z- 
axis, followed by a rotation by —ir/2 around the x-axis. 



The former operation is executed by the phase operator 
P2 = exp(i(7r/2)|2)(2|), while the latter is noting but the 
operation A 2 i. The combination 



= A 21 P 2 U 




(34) 



le 



-1 

27ri/3 
4™/3 



then determines the modulus of n, since 
M|* ) = |0), 

M|* 2 ) = |2>. 

Depending on the outcome |0), |1), or |2) after the mea- 
surement of the qutrit, the number n is divisible by 3 
modulo 0, 1, or 2. The modulus can be found with a se- 
quential measurement scheme: after application of U~ 1 
and measurement of the state |0), the outcome tells that 
the modulus is (if the particle is found in state |0)) or 1 
or 2 (if the particle is not detected) ; in the latter case the 
operator A21P2 is applied and the measurement of well 
|1) will provide the final result, with a modulus 1 if the 
particle is detected in state |1) and a modulus 2 if it is 
not found (then the particle resides in state |2)). 

In the above derivation, we have adopted those phases 
which naturally appear in the simplest manipulation of 
the three- well system — as a result, the matrix Eq. (34) 



is not the canonical (inverse) Fourier transform. The fol- 
lowing steps (for the general case with N counting states) 
then relate the obtained matrix M with the canonical 
form of F _1 . Chosing phases in the definition of eigen- 
states \k) of the counting operator Ci, we interrelate the 
counting and computational bases through 



1 



\k). 



(35) 



Second, let us assume that our physical manipulations 
have produced a measurement operator M such that 



Ml*, 



(36) 



Then the measurement operator M and the canonical in- 
verse Fourier transform F _1 are related via 



with the unitary (phase) operators 



(37) 



/ e tXo 
e lx - 



\ 








(38) 



with phases [x] = [a] , [/?] ■ For the particular situation of 
the qutrit (N = 3), the phases [a] and [0] are defined by 
the operators U p , Eq. (26), and M, Eq. (34). Using these 
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phases in Eq. (37), one easily verifies that the expres- 
sion Eq. ( 34 1 indeed corresponds to the inverse Fourier 



transformation F . 

The generalization of the algorithm to X-qutrit reg- 
isters follows the same steps as above. We assume that 
all elements in the if-qutrit register initially reside in 
the state \0)j, i.e., the state of the register encodes the 
state |0) Q of the computational basis. The qutrits j > 1 
serve the counting of particle clusters: the qutrit j counts 
groups of particles, hence the elementary phase 

shift in Ci is exp(27ri/3 J ). Correspondingly, subsequent 
qutrits are each coupled to the wire a factor 3 less than 
the previous. Let us go through the algorithm: the prepa- 
ration of the j-th qutrit is identical to the first one. Upon 
passage of a particle, the semi-classical states \v) in the 
j-th qutrit pick up phases exp(27rw/3 : '). 

The readout for the divisibility check involves the ap- 
plication of U" 1 to all K qutrits and subsequent test for 
the semi-classical state |0)j in each qutrit: if all qutrits 
j < k + 1 reside in 1 0} ^ and the qutrit j — k + 1 does not, 
then n is divisible by 3 fc . 

In order to recover the number n in the base 3 repre- 
sentation, the qutrit states need to be corrected for the 
passage of incomplete groups of particles before appli- 
cation of M, Eq. (34), except for the first one, which is 
directly measured after application of M. The result 0, 
1, or 2 of the measurement is stored in nx- Before mea- 
suring the next qutrit j — 2, all counters J > 1 need 
to be corrected for the modulus of n: The states \v) 
of the J-th qutrit are given the additional phase shifts 
— 27rnif^l/3 J . Subsequently, the operator M is applied 
to the qutrit j = 2, its state is measured, and the result 
0, 1, or 2 is stored in nx-i- Iterating the process, the 
result 7iK-j+i of the measurement of the j-th qutrit is 
used to correct for the passage of incomplete groups of 3 J 
particles by adding phases — 2-imK-j+iv 3 3_1 /3 J to the 
states \v) of the qutrits J > j. After application of M to 
the qutrit j + 1, its state is measured and the result is 
stored in nx-j- 

Let us briefly analyze the three operators U p , U" 1 , 
and M: since U p merely takes the state |0) into the bal- 
anced state l^o) (but does not act as a Fourier trans- 
formation on the others), its inverse U~ 1 does not de- 
scribe an inverse quantum Fourier transformation. Only 
after augmentation of U~ 1 to the measuring operator M, 
cf. Eq. (34), we arrive at the required inverse quantum 



Fourier transformation allowing for the readout of the 
individual qutrit states. The subsequent readout of the 
qutrit register does not require a fully quantum inverse 
transformation, rather, the semi-classical version using 
sequential measurements and manipulations (executing 
the compensation for the passage of incomplete groups 
of particles) is sufficient. Below, we will encounter other 
implementations, where the preparation operator U p al- 
ready acts as the complete quantum Fourier transforma- 
tion on the semi-classical computational states |n); in 
this case the measurement operator is trivially given by 
M = Ur 1 . The reason for sticking to two different oper- 



ators Up and M _1 in the above discussion is due to the 
simplicity of the preparation step when using U p . This 
is particularly advantageous in the case where one is in- 
terested in the divisibility of n by a power 3 fc , as the 
readout only involves the inverse operator U~ . Using 
the full Fourier transformation M _1 for the preparation 
instead leads to a much more difficult hardware imple- 
mentation of the preparation, cf. the following section 

E3 



A. Generalization to qudits 

The further generalization to base-d counting with qu- 
dits follows the same ideas as those developed for the 
base-2 and base-3 counting with qubits and qutrits. In 
order to set up the algorithm, we have to define the three 
steps 'preparation through U p ', 'counting with Ci', and 
'measurement' with the inverse Fourier transformation 
M. The first two steps are clear: the preparation of the 
initial balanced counting state l^o) starts from the com- 
putational state |0) and proceeds with the subsequent 
shift of weight in the wave function to the neighboring 
well, always leaving behind an amplitude with weight 
1/Vd; the individual steps involve the operators 



(I 

... 
... 

\o 



\ 



cos x 1 sin x 
zsinx cosx 



(39) 



I) 



with the (k — l,k) non- trivial 2x2 block shifting the 
amplitude between the wells k — 1 and k, cf. Eq. (20). 



This operation is physically implemented through lower- 
ing the barrier between two neighboring wells k — 1 and 
k. Similarly, the counting operator Ci is given by the 
straightforward generalization of Eq. (27) involving the 
rotation- or phase-operators 



/I 



1 

n <M> 



\ 

... 
... 



\ 



(40) 



lj 



with <fi = 2nk/d and implemented through changing the 
gating bias of the well k. 

The implementation of the inverse Fourier transforma- 
tion M with the help of the physical operators Eqs. (39) 



and (40) is more difficult but still possible, since this set 



of operations (gates) is universal, i.e., given a unitary dxd 
matrix, it can be constructed from a product of operators 

The proof of this statement 



made from A£ , 



-i ^d P£. 
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is similar to the proof of the universality of two-level uni- 
tary gates, cf. Ref. 9: Given an unitary d x d matrix 
(or operator) U, the idea is to reduce U in an iterative 
procedure to the unit operator by right-multiplication 
with amplitude and phase operators. E.g., consider the 



entries U 



o.fc-i 



= P = 



■'■■pp 



and Uo.fe = a = \a\e l 



in U (we number rows and columns with indices from 
to d — 1). The product U Pf_i A? fc l generates the 

new entries j3 e 1 ^ cos \ + ia sin \ and a cos \ + i/3 e 1 ^ sin x 
and we can replace the new entry at the position (0, k) 
by if we choose the angles x = arctan(|a//3|) and 
4> = tt/2 + <p a - <pp (if /? = then x = 7T /2, = 
and a — implies zero angles). The new (0, k — 1) entry 
reads j3' = ie ltPa \f\a\ 2 + |/3| 2 . Repeating this step d — 1 
times we replace the top row by zeros, except for the (0, 0) 
entry, which we bring to unity with an additional phase 
operation. Hence 2c? — 1 elementary amplitude and phase 
operations take the first row to the vector (1, 0, 0, . . . , 0) 
and a total of — 1) = d 2 elementary operations 

take the unitary operator U to unity; the desired opera- 
tor then is obtained by a simple inversion. 

In order to illustrate the procedure, we derive the 
Fourier transform M _1 , 



M" 




a 2m/3 



1 

4tti/3 



47ri/3 „87ri/3 



(41) 



for the qutrit. In the implementation of the first step we 
find the angles x — arctan(l) = 7r/4 and <p = tt/2 and 
the product M x = M~ P ± A 21 produces the matrix 



, . iV2 \ 

Mr 1 = -7= 1 -VV2 -iy/3/2 
V6 ^ 1 -i/V2 J 



(42) 




The angles for the second step read x — arctan(v / 2) and 
(j) = tt and we obtain the matrix M^ 1 = M^ 1 PJ A^ PJ 



(43) 



Next, the angles for the third step read \ — tt/4 and 
4> = tt/2 and the product M3 1 = M^ 1 P^ /2 A^ 4 gen- 
erates the second row in the form (0,1,0) already — no 
further phase rotation is needed in this step. Finally, 
we have to compensate for the phase of the (2, 2) ma- 
trix clement, which is done by the additional rotation 



M 3 ip;^/ 2 1 and we arrive at the unit matrix. Col- 



lecting all factors, we obtain the measurement operator 
M expressed through elementary shift and phase opera- 
tors, 



= K 



tt/2 .tt/41 
M 2.1 J 



pTT ax p^irp 71 
r o rt i,o r oJL r i 



77/2 .7T/4 



A 2,l ] [ P 



37r/2i 
2 \ 



(44) 



The result Eq. (44) corresponds to the expression Eq. 



minimal algorithm with a preparation step U p generat- 
ing the first computational state with phases (1, i, — 1) 
instead of the canonical ones [the trivial phases (1, 1, 1)]. 
As a result, the final measurement operator Eq. (34) cor- 



(34) up to phases. Indeed, above, we have presented a 



responds to the inverse Fourier transform up to phases. 
Of course one could easily introduce additional phase op- 
erations P^ and remove the non-canonical phases (with 
Up = P^/ 2 P£ A^ with x = arctan%/2), however, 
such additional gates only render the algorithm more in- 
volved. Also, it is important to note that the divisibility 
check does not require the implementation of the inverse 
Fourier transformation and hence it is worth while to 
know how to implement a minimal preparation operator 
U„. 



VI. IMPLEMENTATIONS OF QUTRITS AND 
QUDITS 

The base d counting and factorization algorithm obvi- 
ously requires a set of suitable and well operating elemen- 
tary quantum devices, qubits, qutrits, or qudits. Starting 
out from qubits and their analogy with a spin-1/2 system, 
the most natural generalization is to try a spin-1 system 
for the implementation of qutrits and possibly a spin-<i 
system for the qudits; although this idea can be realized 
in principle for the case of a spin-1 system, the prepara- 
tion and measurement algorithm is rather complex (see 
below) and should be viewed as a Gedanken experiment 
rather than a realistic proposal. The next idea then is to 
generalize the concept of the double-dot charge qubit — 
this road has been pursued above and works fine in the- 
ory, however, the implementation of multi-dot charge 
qubits may turn out difficult. As an alternative, one may 
try to emulate a d-spin qudit (<i-lcvcl system) by a system 
of spin-1/2 qubits (two-level systems). Provided we ad- 
mit two-qubit interactions in our manipulation scheme, 
we then are able to define all the necessary operations 
required by the algorithm. Hence, we can offer a scal- 
able route for the implementation of qudits using qubits 
as elementary units. Below, we will discuss the various 
issues related to the implementation of qutrits and qu- 
dits in more detail, putting our main emphasis on the 
understanding of the qutrits, their implementation and 
manipulation. 



A. Spin-1 qutrit 

The most straightforward attempt to generalize the 
base 2 counting with qubits (spin-1/2 two-level systems) 
to the base 3 counting with qutrits makes use of a 
spin-1 three-level system with the orthogonal (compu- 
tational) basis \l,m) z — |l,l) z = |0), \l,0) z = |1), and 
|1, — 1) 2 = 1 2), where I and m denote the angular momen- 
tum and magnetic quantum numbers. As in the previous 
cases (two-level qubit, three-level qutrit), we have to pre- 
pare the system in the initial counting state e.g., the 
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lowest harmonic of the computational basis or another 
balanced state. The logic taking us to a valid balanced 
state is the following: i) we first note, that experimen- 
tally we can produce either axial states (polarized along 
a direction n) such as |1, l) n , or planar polarized states 
1 1, 0) n , hence we will focus on these types of states, ii) As 
one cannot construct axial states with equal weights for 
all basis states |0), |1), |2), we concentrate on the planar 
states. Hi) the most general planar state (with a direc- 
tor parallel to n defined by the direction angles tp and 9) 
assumes the form (in the computational basis) 



z 


1 






/ 


1 ^0 





















|l,0) n 




-smt/e ' v 
y/2 cos 9 



(45) 



We demand that all components have equal weights, 
hence 9 = arctan-\/2; choosing ip = 7r/4, we obtain the 
balanced state 



-i7r/4| 



ID 



|2) 



(46) 



the planar state |l,0) ec with a director e c = (1, 1, l)/v3 
pointing along the body diagonal; other values of tp corre- 
spond to another direction e c rotated around the z-axis. 

The other counting states then are obtained by rotat- 
ing |^o) by the angle 27r/3 (anti-clockwise) around the 
z-axis, = U z (2n/3)\* ) and |tf 2 ) = E/,(27r/3)|*i); 
physically, this rotation is achieved by ensuring that pass- 
ing electrons create a local magnetic field pulse along the 
z-axis. To simplify our further discussion, we exchange 
the two axes, the one for the preparation and the one 
defining the counting field: our task then is to generate 
an initial planar polarized state |1, 0) 2 and implement the 
counting step through rotation by an angle 27r/3 (anti- 
clockwise) around the axis e c . 

While the creation of axially polarized states is rather 
straightforward, creating a planar polarized state |l,0) n 
is non-trivial. Indeed, while we can make use of a simple 
deterministic procedure to prepare an axially polarized 
state (by switching on a magnetic field and relaxing the 
spin through coupling to a bath), the planar polarized 
state required here is more difficult to obtain. A simple 
preparation can be implemented with the help of a Stern- 
Gerlach apparatus, however, this procedure is a statisti- 
cal one, with a one-half probability for a positive out- 
come. In fact, preparing an initial state polarized along 
x, 11,1), = + \l,-l) z )/2+ |1,0) Z /V2, using a 

Stern- Gerlach setup directed along z, and selecting par- 
ticles with a straight trajectory (i.e., selecting the 'middle 
spot' in the Stern-Gerlach setup) we obtain the desired 
spin state |1, 0) 2 (a procedure to generate a proper initial 
state with unit probability is discussed later). 

With the counting state |^o) = 1 1 , 0} ^ properly pre- 
pared, we define the counting step through a rotation by 
the angle 27r/3 around the axis e c = (1,1, l)/v3, cf. Fig. 



FIG. 7: Spin-1 qutrit. After preparation in the state |1,0} Z , 
the system is in a state of planar (xy) polarization, i.e., the 
spin has no component along z. Upon passage of a particle, 
the polarization plane is rotated (by the angle 2n/3 around 
e c ) from xy to yz, for the second particle to zx, and for the 
third back to xy. The divisibility test involves a measurement 
with a Stern-Gerlach apparatus with the field directed along 
the z-axis. 
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Ci = cxp[— i(27r/3)e c • "L/h] 

1 e-"/ 4 



(47) 



exp 



^= I e"/ 4 e~™ 
3v3 \ q e i7r/4 _i 



*0> = 


|1,0)„ 


*1> = 


|1,Q>«, 


*2> = 


11,0),. 



As the three-fold iteration takes us back to the original 
state |1, 0) Zl our counting operator is cyclic and generates 
the complete orthogonal counting basis 



(48) 



The measurement step for the divisibility check af- 
ter the passage of the particles involves a second Stcrn- 
Gcrlach experiment directed along the z-axis — if the par- 
ticle moves again on the straight trajectory, its polar- 
ization was unchanged by the passage of the n parti- 
cles and hence n is divisible by 3. The other states 
11,0)* = (|W),-|MU)/>/2and |1 > 0>„ = *(|1,1>» + 
cf. Eq. (45), give no contribution to the 



|1,-1>,)/V2, 

signal on the straight trajectory, cf. Fig. [HJ 

As usual, the measurement of the number's modu- 
lus (the counting measurement) is more involved. We 
then make use of the other outcomes |1, l) z and |1, — l) z 
of the second Stern-Gerlach experiment and bring them 
to interference further down their trajectories, cf. Fig. 
[8] Testing the resulting state, e.g., the state |l,0) y = 
j(|l,l) 2 + |1, — l) z )/v / 2 (we choose symmetric trajecto- 
ries with equal phases) in a third Stern-Gerlach appara- 
tus polarized along the y-axis then identifies the counting 
state 1 1, 0) y on the straight trajectory (the middle spot, 
the number's modulus is 2); if no spin is measured, the 
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counting state |1, 0)^ has been realized and the number's 
modulus is 1. 



FIG. 8: Setup for the counting with a spin-1 qutrit. Send- 
ing the initial polarized state |1, 1}^ into a first Stern-Gerlach 
(SG) apparatus with axis z and selecting the middle (un- 
deflected) trajectory, we obtain the planar state |1,0} Z with 
probability 1/2. The counted particles passing by rotate this 
state around the axis e c = (1, 1, l)/\/3 and generate the pla- 
nar counting states |1,0} Z , 1 1, 0}^ , and |1,0} H (note that the 
counter particle has to be stopped and trapped during the 
time of interaction with the particles passing by). These are 
analyzed in a second SG apparatus with axis z — if the spin 
is not deflected, the number n of counted particles is divisi- 
ble by 3. Combining the two deflected beams and analyzing 
the superposition in a further SG apparatus directed along 
the y-axis, we find the mudulus of n, which is 2 if the spin is 
undefiected and 1 else. 

Finally, we look for a preparation step with unit effi- 
ciency. This can be achieved by making further use of the 
states |l,±l) z behind the first Stern-Gerlach device (we 
still use the input state |1, l) x to the Stern-Gerlach ap- 
paratus): superimposing these two states, we obtain the 
component l) z + 11, — 1) 2 )/2 = |1, Q) y /^/2 and rotat- 
ing this planar state back to the z-axis, we have managed 
to transform the initial spin |1, l) x into two subsequent 
wave-packets with weight 1/2, each describing the same 
spin state |1,0) 2 . Hence, although we cannot generate 
a planar state out of an axial state by simple rotation, 
the use of a Stern Gerlach apparatus and proper manip- 
ulation of all three amplitudes of the split wave function 
allows one to generate the desired planar state, though 
split into two wave-packets with weight 1/2 each. 

Obviously, the above scheme is a very complex one and 
should be regarded as a 'Gedanken' experiment rather 
than a realistic setup. Nevertheless, it is interesting to 
see, that a spin-1 qubit can, at least in principle, be used 
for the implementation of the base-three counting algo- 
rithm. 



B. Triple-dot qutrit 

The triple-dot qutrit has been discussed above in sec- 
tion |Vj Preparation, counting, and readout can be prop- 
erly implemented via voltage pulses acting on the semi- 
classical states (phase shifts) or on the barriers in be- 
tween (amplitude shifts). A drawback is the need to 
design a new device when going from base 2 to base 3 
counting and, more generally, each time a new prime fac- 
tor is to be tested. The emulation of qutrits (and qudits) 
through qubits described below allows one to stay with 



only one computational unit for all base-d counting tasks 
and factorizations. 



C. Emulation of spin-1 qutrit 

The basic idea we pursue here is to emulate the qudits 
needed in the base-d counting and factorization through 
simpler qubits. We start by combining 2 qubits into a 
qutrit. 



1. Emulation using a spin triplet 

An obvious way to choose three appropriate states in 
the product Hilbert space 'Hi/2 ( X ) '%i/2 of the two qubits is 
to make use of the decomposition into singlet and triplet 
sectors, 



H 1/ 2®H 1/2 =H (BH 1 , 



(49) 



and use the three-dimensional triplet space Hi. Con- 
trary to the simple spin-1 qutrit, the emulated version 
using two qubits provides us with the necessary degrees 
of freedom to perform all of the required steps (prepara- 
tion with Up and readout with M) in the counting and 

factorization algorithm with the help of one- 
and two-qubit (a^a^) operations. Referring to the pre- 
vious paragraph, we start from the computational basis 
|0> = | tt>, |1) - (I U) + Ut»/V2, and |2) = | W) 
and seek for those manipulations which provide us with 
the counting states in the form of three orthogonal pla- 
nar polarized states with balanced weights. The latter 
take the form Eq. (45), whereas the axial states can be 



written as (again in the computational basis) 



IU>. 



/ (1 + cos 6) e-tf 

V2 sin9 
\ (l-cos6»)e^ 



(50) 



with n the direction of the spin axis. Note that once 
these states are chosen, the operators <Jz and a l ^ ) crf ) 
only change the relative phases in |l,0) n and |1, l) n but 
not the relative weights in the amplitudes. It is then our 
task to construct three orthogonal planar states out of 
the computational basis, which consists of two axial and 
one planar state. In this way, we can find a preparation 
operator U p which is congruent to the Fourier transfor- 
mation, i.e., Up = M _1 transforms the computational 
basis {|A:)}|^ =0 into the counting basis {l^n)}!^)- 

We first note, that the planar states Eq. p5 1 always 
have the same weight in the components |0) and |2) and 
hence we first have to rotate the axial states of the com- 
putational bases into the xy-plane. This is done with the 
operator e l7ro W 2 where a x — (er!-. 1 ' + a^ } )/2 and results 
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(up to a phase) in the states [9 = ir/2 in Eq. (50) 



1 



1 



|l.l) e „ = - | V2^ 



1 



-1 



|1. I.) . = 2 [ ^< |. (51) 



cf. Fig. [9j Note that the simple rotation takes the planar 
state |1) into a planar state (with a director along y). In 
order to map the axial states in Eq. ( |5"Tj | to planar states 
we choose 6 = n/4 in Eq. ( [45] ) and obtain the candidate 
states 




(b) 



y 



X 



11,0), 



V2 



(52) 





!*2> 








f/ % 




v 







the planar and axial states in Eqs. (52 1 and (51) now 




have equal amplitudes. In order to select the appropriate 



phases ip in the planar states Eqs. (52 ), we note that the 
phase differences between the |0) and |2) components of 
the axial states are equal to ±7r; these have to be matched 
with the phase differences tt — 2ip in the planar states and 
hence we choose planar states with ip = 0, 7r, 



|1,0) 3 




|l,0)x,* = 




(53) 



The two states are characterized by directors pointing 
along the irz-diagonals, cf. Fig. [9j In order to map the 
axial states Eq. (51) to the planar states Eq. (53) we 



X 



FIG. 9: Rotations generating the quantum Fourier transfor- 
mation on the triplet sector of a two-qubit system; we denote 
the axial states 0) — ft) and 2) = by double arrows 
and planar states (e.g., 1) = [ f4-) + I by a director, 

(a)— >(b) Rotation by —n/2 around the x-axis, (b)— »(c) con- 
ditional rotation transforming the two axial states into pla- 
nar ones with directors along the diagonals in the £z-plane, 
(c)— >(d) rotation by <p = — arctan(l/\/2) around the z-axis 
of three planar states to position them symmetrically around 
the z-axis. The last step defines the counting states which 
transform among one another through a rotation by the an- 
gle — 27r/3 around the z-axis. 



make use of the two-qubit operator 



(planar) states \^ n ) = U p |n), 



lb 



(54) 



l*o> 
l*i) 



3 -5iri/12 
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One can confirm, that the counting step involving the 
application of the operator Ci = exp(27ri cr z /3) with 
(j z = (a z 1} + a { z 2> )/2, indeed transforms the three counting 



The latter adds the phases x i~X) to the components |0), 
1 2) (|1)) and hence leaves the relative phase between the 
components |0), |2) unchanged while adding a relative 
shift — 2x to the middle one. Hence choosing \ = — 7r /4, 
we can map the axial state |l,l) e to the planar state 
IDO)^ and |l,l)_ e to the planar state |1,0) X) _ Z , cf. 
Fig. [9j Since the component |1) in the planar state 
|l,0)e y has weight 0, the operator leaves it un- 

changed. 

Hence we have arrived at three orthogonal planar 
states directed along y and along the xz-diagonals; a 
final rotation e 1 ^* with a x = (a x 1} + er<, 2) )/2 by <p = the following paragraph, section |VI C 2\ 
arctan(l/v2) around the x-axis then arranges the three 
states symmetrically around the z-axis. Combining the 
three operations, we obtain the preparation and inverse 
measurement operator 



-|2>, 



e -i57r/12 



(56) 



|2). 



states in Eq. (|56|) into each other, i.e., |* n ) = C™|4'o) 



Cn|*o) (note the symmetric definition of a z used here 
and compare with the corresponding definition of Ci in 



2. Efficient emulation 



u, 



vp<T x 



-inert 1) o-W /4 i va - x /2 



(55) 



The actions of the three operations in Eq. ( 55 1 on the 
computational basis states {|^)}|^ = o are illustrated in 
Fig. [9j the resulting counting basis is given by the three 



Although the emulation of the qutrit by the triplet sec- 
tor of two qubits can be done consistently, this scheme is 
inefficient in scaling up to larger primes. Indeed, emulat- 
ing a spin-2 system within this scheme requires 4 qubits, 
whereby most of the Hilbert space (an eleven dimensional 
sector) is not used since only the component %2 out of 
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the decomposition 



H i/2 - H z 



(57) 



is needed. A suitably scalable case should make maxi- 
mal use of the emulating qubits, i.e., use a large fraction 
of the Hilbert space. In choosing the appropriate sub- 
set of states in the multi-qubit Hilbert space, we have 
to select states with an equidistant spectrum. E.g., for 
the case d = 3 and 2 qubits, we choose three consecu- 
tive states from |0> = | ft), |1) = I t4->, |2) = | It), 
1 3) = | ||) and an elementary phase shift operator 
Ci = cxp[—i\(a l z 1) + <r ( / } /2)]; upon passage of a parti- 
cle, the four states then pick up the phases exp(— 3iA/2), 
cxp(— i\/2), exp(iA/2), and exp(3iA/2). Choosing A = 
27r/3 and the first three states |0), |1), and |2), the (rel- 
evant sector of the) operator Ci, 
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(58) 



coincides with the expression Eq. (|27j), up to an overall 
phase exp(3iA/2) = exp(i7r) [to be added to ( 58 ) ] . The 



generalization of this scheme to an _E-qubit emulation of 
a qudit is straightforward, 



exp 



(i-i) 



i=i 



(59) 



In addition to the shift operator C\, we need to know 
the form of the preparation operator U p . For the above 
qutrit emulation, this operator takes the form (see Ap- 
pendix [A] for the derivation) 



with the angle 9 = 2 arctan(l/-\/2). As this is not yet the 
full Fourier transformation on the computational basis 
\k), k = 0, 1, 2, we still have to find the measurement op- 
erator M . In order to accomplish this task, one notes that 
the three counting states |\& n ), n = 0, 1, 2 are entangled, 
while the three computational states \k), k = 0, 1,2 are 
not. The inverse Fourier transformation M then has to 
disentangle the counting states, a criterion which helps 
us in finding its explicit form. In the end, the measure- 
ment operator M is obtained in the form of a product of 
three unitary operations Uj, i = 0, 1, 2, 



= UoUiUs, 



(61) 



where U2 and Ui serve to disentangle the three counting 
states |*n), n = 0, 1, 2, and Uo is a conditional Hadamard 
operation, turning the spin of qubit 2 into the z-axis if 
the qubit 1 is in the state | f). The detailed derivation 
of Eqs. (60 1 and ( |6T] ) and the form of the operators Uj, 
i = 0, 1, 2, is given in Appendix |A"| 



VII. RELATION TO PHASE ESTIMATION 
ALGORITHM AND ITS APPLICATION 

It turns out, that our counting algorithm has much 
in common with the phase estimation algorithm (PEA); 
the following discussion of the PEA is formulated in a 
way as to make this connection apparent. The phase 
estimation algorithm first appeared as a part of Shor's 
factorization algorithm^; an extended separate algorithm 
was presented by Kitae\W and later by Cleve et alP^. 
The phase estimation algorithm attempts to find the 
'phase' < <p < 1 in the eigenvalue exp(2mip) of a 
unitary operator U associated with a given eigenvector 
\u). In the version of Refs. I9|12|, this is achieved with 
the help of two qubit-registers, one of which (the sec- 
ond) is storing the vector \u) and acts on it with the 
operators U 2J \ j = l,...,K, to generate the phases 
exp(27ri 2^ 1 p). The other (first) register consists of K 
qubits and produces the desired phase estimate in the 
following manner: with all qubits initialized in the state 
|0) j , a Hadamard operation generates the balanced states 
(|0)j + |l)j)/V2 for all qubits in the first register. A con- 
trolled U 2 1 operation between the second register and 
the j-th qubit in the first register then puts the qubit 
into the state (|0)j+exp(27ri2 A_: '(^)|l) :) )/v / 2, hence gen- 
crating the quantum Fourier transform (we express the 
product state through the computational basis and as- 
sume that the phase (p can be represented by K binary 
digits) 



F(|2*¥>» = 
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(62) 



in the first register. A final inverse Fourier transforma- 
tion then generates the state \2 K ip) Q \u) and the projec- 
tive measurement of the if-qubit register in the com- 
putational basis provides us with the phase ip; for an 
arbitrary phase < ip < 1 we obtain a X-binary-digit 
estimate t^dig of the phase ip. 

Comparing this algorithm with our counting setup, we 
identify the action of the n particles traversing the quan- 
tum wire with the action of the second register in the 
PEA, with the correspondence |n) s ■<-> \u). The con- 
trolled U 23 operation in the PEA is replaced by the cou- 
pling of the wire to the qubits: the interaction of the n 
particles with the last (K-th) qubit has to be identified 
with the action of the controlled U operator in the PEA, 
hence ip = n/2 K . The qubits j < K are more strongly 
coupled to the wire, that corresponds to higher powers of 
the operator U in the PEA; in fact, the j'-th qubit cou- 
pling is enhanced by the factor 2 K ~i and its interaction 
with the particles in the wire corresponds to the action of 
the controlled U 2 3 operator. Finally, the intermediate 
states in Eqs. Q and ( [62] ) agree with one another with 
the identification n/N = n/2 K -H> cp. The final states ex- 
hibit the correspondence |n) Q |n) # «-> \2 K Pdig) q\u) . Note 
that our divisibility algorithm has no counterpart in the 
PEA. 
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This analogy imm ediat ely allows us to profit from the 
performance analysis^^ of the PEA: Assume that we 
wish to measure the phase <p in the PE problem to an 
accuracy of 1/2 A (i.e., we want to encode ip with A bits) 
and be sure of our measurement result with a probability 
P = 1 — e at least, then the setup must involve K — 
A + [log 2 (2 + l/2e)] qubits. 

This result can be applied to our counting algorithm. 
Consider the case, where a non-integer number x = n+Sn 
has passed the counter, with n an integer and < Sn < 1 
a real number. Such a situation may occur when the 
interaction between the particles and the counting qubits 
is still finite at the moment when the readout procedure 
starts, corresponding to the passage of a fraction of a 
full charge. Concentrating first on the implementation 
with a fully-quantum inverse Fourier transformation, the 
performance analysis of the PEA tells us that we still can 
measure the number to any desired precision. E.g., if we 
want to be able to measure the number n < N — 2 K such 
that |n meas — x\ < 1/2 with a probability P = 1 — 2~ r , 
we need to be able to resolve fractional charges Sn ~ 
2 _r <C 1, i.e., we have to add additional qubits which 
measure half-charges (turning by 2ir on the passage of 
one particle), quarter charges (rotating by 4tv), etc. The 
entire setup then has to involve w if+log 2 (l/2~ r ) = K+ 
r qubits. This result can be extended to qudits: requiring 
a precision P = 1 — d~ r , we need rj K + log d (l/d _r ) = 
K+r qudits. Hence, we can trade additional qudits in the 
counting process against a higher probability to obtain a 
correct integer result. 

Next, we consider the semi-classical inverse Fourier 
transformation; we demonstrate below that this semi- 
classical scheme exhibits the same stability as the fully- 
quantum version, although the passage of a fractional 
charge strongly affects the conditional measurement of 
subsequent qubits, e.g., for Sn = 1/2 the measurement 
of the first qubit gives a random input for the direction 
of measurement of the second qubit. However, this er- 
ror does not propagate through the entire measurement 
scheme. Instead, the measurements of higher qubits re- 
cover from false results measured for lower qubits. For- 
mally, this can be proven by comparing the two prob- 
abilities Pqp(n]x) (using a full quantum Fourier trans- 
formation) and P sc p(n;x) (using a semi-classical Fourier 
transformation) to find the integer number n when a non- 
integer number x has passed by the qubit register. 
Indeed, let us calculate the probability P bc f(ji, x) to 



find the integer result 



K 



upon mea- 



surement of a non-integer signal x by the qubit regis- 
ter. The probability pi(nj(,x) that the first qubit pro- 
vides the value Uk is given by the matrix element be- 
tween the qubit state (|0)i+e 27 " x/2 |l)i)/%/2 and the state 
(|0)i + e 7rmx \l)\)/\/2 to be measured (i.e., the projection 
| + y)x for nj{ = or | — y)\ for ur = 1, cf. Sec. Ill) 



Pi(n K ;x) 



(l|)(|0) 1 +e 2 ^/2|i )l) |i 



Next, the conditional probability to measure the value 
riK-i for the second qubit is given by the product 
p 1 (n K ]x)p 2 {n K _ l ,n K ;x) with 

p 2 (n K -i,n K ;x) = |( 2 (0| + e -*«"*-i+W2) 2 <i|) 

(|0) 2 + e 2 ™/ 4 |l) 2 )| 2 /4. (64) 

Using Eq. ([8]), we can rewrite the first factor in 
pi(jiK-\, nx', x) in the simpler (and equivalent to the sec- 
ond factor) form 

. |( 2 (0|+e- 2 ™/4 2(1 | )( | ) 2+e 2— /4| 1)2 )|2 
P2{n]x) = . 

(65) 

The straightforward iteration of this scheme then pro- 
duces the final result for P sc -p(n,x) in the product form 



A" 



P scF (n,x) = Y[pj(n;x), 

3=1 



(66) 



where 



Pj(n,x) 



| ( . (0| + e -2 7 rW^. ( i| )( |o) J . + e 2^| 1) . ) p 



(67) 

Evaluating the product, the result Eq. (66) is easily 
rewritten in the form P sc p(n,x) = \q^ti\W x )q\ 2 with 
the generalized counting state (cf. Eq. (pi)) 



(63) 



2"-l 

|* x ) = -= exp(2inxk/2 K )\k) Q . (68) 
^ 2 fc=0 

Expressing the counting state I^^q as a Fourier trans- 
form of the computational state |n) Q ), | x I , „)q = F(|n) Q ), 
we arrive at the result P sc p(n,x) = Iq^I^qI 2 , where 
\x) Q = F _1 (| v E , 2;)q) is defined in terms of the back trans- 
formed counting state. But this is nothing else than the 
probability Pqp(n,x) to find the number n in a one-shot 
measurement of the qubit register after application of an 
inverse quantum Fourier transformation on the detected 
state \^ x )q, hence P SC F(n 7 x) — P q F(n,x). We conclude 
that the semi-classical and the fully-quantum Fourier 
transformation exhibit the same stability to systematic 
errors introduced by incomplete (non-integer) counting. 
Furthermore, we note that the semi-classical algorithm 
is rather robust with respect to random errors; the latter 
can be handled with a classical multi-qubit error correc- 
tion scheme combined with a simple majority rule, cf. 
Ref . H 



A. Quantum metrology: voltage measurement 

The above discussion allows us to use our counting de- 
vice to measure continuous variables. The insight paves 
the way for its use as a quantum voltage-detector, a par- 
ticular form of an analog to digital converter (ADC). 
Consider a setup similar to the one in Fig. [T] but with 
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FIG. 10: Quantum voltage-detector made from K counter 
qubits. A voltage V applied over a time interval t generates a 
continuous phase <p = eVt/h, which can be measured by the 
iGqubit register and translated into a binary digital output 
signal. 



the wire replaced by a finite metallic object, see Fig. 10 



Assuming that the applied voltage translates into a 
relative shift 8k = axeV (<5i = a\eV) of the most weakly 
(strongly) coupled qubit [which is the qubit j — K (j — 
1)], the reading of the qubit array after detaching the 
voltage provides us with a binary digital number n which 
translates into the phase cp according to 



tp ■ 



eVt 

IT 



27r n 



2ir n 



(69) 



provided that the quantity 2ot\eVt/h can be represented 
by an integer number. In order to handle successfully the 
general case with an arbitrary drive Vt, we have to add 
additional qubits in order to arrive at a digital estimate of 
ip: for a result with a relative precision 1/2 and with a 
probability better than 1 — e, we need a measuring device 
with K + [log 2 (2 + l/2e)] qubits. 

The sensitivity of our quantum ADC improves as 1/t, 
where t is the time of observation; this is because the 
accumulated phase increases linearly in time, while the 
most sensitive qubit always resolves phases of 7r. The 
precision then is by a factor \j\ft better than the usual 
classical sensitivity, which scales as and agrees 

with standard expectations^. In particular, using the 
straightforward algorithm described in Sec.|ll]to measure 
a voltage V < V max = irh/er with an accuracy 8V, we 
either need M = (7r/</>) 2 = (V max /SV) 2 qubits or have to 
run repeated experiments over a time t — Mt, where r is 
the measuring time of the individual experiment. Solv- 
ing for the desired accuracy 8V , we obtain the scaling 
8V ~ y/hV max /et, i.e., the precision improves only with 
the square-root of the overall measuring time t, the same 
as for the usual classical case. 



VIII. MULTI-PARTICLE ENTANGLEMENT 

Another application of our counting device is the gen- 
eration of multi-particle entangled states with the help 
of a Mach-Zehnder interferometer, see Fig. [TT] Inject- 
ing particles into the device through the lower left arm, 



the splitter generates a superposition of number states 
in the two arms of the interferometer; measuring the 
counter placed near the upper arm and selecting a par- 
ticular reading projects the system to the desired entan- 
gled state. The functionality of the device has been de- 
scribed in detailed in Ref. [4] before. Here, we use the 
device with our qutrit counter to generate the original 
GHZ (Greenberger-Horne-Zeilinger) state^ and to un- 
veil in more detail the entanglement between the counter 
states and the physical number states in the quantum 
wire; the generalization to other cases with more parti- 
cles and counters follows the previous discussion in Ref. 

H 

Consider a particle entering the Mach-Zehnder inter- 
ferometer from the lower-left lead and propagating along 



one of the two leads U or D, see Fig. 11 The wave func- 
tion can propagate along two trajectories, the upper arm 
U where the particle picks up a phase <pu and the counter 
is activated, or the lower arm D accumulating a phase ipo 
and leaving the counter state unchanged. The total wave 
function evaluated at the position A then assumes the 
form 

*ia = te l ^\ ft) ® |*i)+re i,w> |^) ® |* ), (70) 

where t and r denote the transmission and reflection co- 
efficients of the beam splitter and we have introduced a 
pseudo-spin notation to describe the propagation of the 
particles along the two arms: a pseudo-spin ft (JJ.) refers 
to the particle propagating in the upper (lower) arm. The 
qutrit state depends on the particle's trajectory and reads 
either \^!\) if the particle has passed in the nearby upper 
arm or |^q) if t ne particle passed through the lower arm 
of the interferometer. 




FIG. 11: Mach-Zehnder interferometer with qutrit counter. 
Particles enter the interferometer through the left leads (here 
the bottom lead) and are measured on the right. The qutrit 
counter in the upper arm U detects the passage of particles. 
The magnetic flux "!> through the loop allows to tune the phase 
difference when propagating along different arms. 



Next, we inject three particles (from the bottom left) 
into the Mach-Zehnder (MZ) loop. We assume the three 
individual wave functions describing the initial state to 
be well separated in space, allowing us to ignore exchange 
effects in our (MZ) geometry. The wave function at the 
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position A then reads 



+tr 2 e ^u+2 VD )^ AA) + ^ AA) 

%{2vu+tpD) [\1\AA) 



^^^'[niD + itit) 



l*o) 
(71) 



Assuming scattering coefficients for a symmetric beam 
splitter, e.g., tt* = 1/2, r 2 = (-l)/2, and rt* = 
±ie~ tipt /2 (with ipt the transmission phase), the pro- 
jection to the counter state |^o) provides one with the 
GHZ-like state |* GHZ ) = {^A^Tie^^-^-^MA 
, JJ-))/v2; manipulation of the flux $ in the Mach-Zehnder 
loop then allows one to implement the desired entangled 
state. Furthermore, the wave function Eq. (71 1 unveils 



the entanglement between the (split) number states and 
the counter. Note that the indistinguishability of par- 
ticles exploited in the above entanglement process is an 
'artificial' one defined by the qutrit detector, rather than 
the 'fundamental' one of identical particles. 



IX. DISCUSSION AND CONCLUSION 

Summarizing, we have generalized the binary (base 2) 
quantum counting algorithm and divisibility test by 2 k 
to ternary (base 3) and higher counting systems. This 
extension is quite non-trivial in several respects: On the 
device level, the qubits used in the base 2 algorithm have 
to be replaced by qutrits for a ternary counting system 
and to qudits for a base d algorithm. Since the algorithm 
is based on two subsequent quantum Fourier transforma- 
tions, suitable manipulation schemes have to be defined 
in order to implement a quantum Fourier transformation 
on the level of individual qubits, qutrits, and qudits. Fur- 
thermore, rather than developing new hardware for every 
new counting base, we have discussed how to use qubits 
in order to emulate qudits, with particular emphasis on 
the qutrits. 

Starting from binary counting with qubits, it seems 
not immediately clear how to generalize the concept. It 
turns out, that defining the quantum counting task on 
an elementary level through a one to one correspondence 
between the counting objects and distinguishable states 
in a Hilbert space provides us with a constructive scheme 
how this task can be achieved. Also, the analysis of the 
unary counting scheme naturally introduces the quantum 
Fourier transformation as the basic operation in a non- 
demolitian counting process. Indeed, the task of count- 
ing naturally introduces a shift operator Ci in the count- 
ing space, taking one counting state I'ffjq into the next 
\^n+i)<3- When expressing these counting states through 
the eigenstates \n) Q of Ci, then the counting operation 
only adds a phase exp(27rm/7V) to each of these states. 
Hence using these eigenstates \ri) Q , which are nothing 
but the Fourier transforms of the counting states \^ n ) Q , 



as our computational basis provides for us a 'soft' non- 
demolitian counting scheme. Choosing another compu- 
tational basis involves an energy exchange between the 
counted object and the counting system and introduces 
a much more severe perturbation. 

The aforementioned basic understanding of quantum 
counting has provided us with a constructive scheme for 
the counting algorithm: starting out with a set of mea- 
sureable states (the computational basis {|«)q}) which 
evolve with a prescribed phase accumulation upon the 
passage of particles, we have to prepare out of them a 
balanced state \^o)q which is used as the first counting 
state. With the appropriate phase increment picked up 
during counting, this first counting state evolves to the 
next and returns to the first one after an iV-cycle. As we 
have seen, we do not have to enter the cycle in the low- 
est harmonic, any harmonic will do, and even arbitrary 
phases in a balanced state (an equal-weight superposi- 
tion of the computational basis states) are acceptable. 
In the latter case, the Fourier transformation is modi- 
fied with additional phases, which are not harmful to the 
algorithm, however. 

Equipped with this general scheme, we have been able 
to generalize the binary algorithm to a ternary and 
to base-ci counting. We have seen, that one possible 
straightforward extension of the hardware, going from 
a spin-1/2 to a spin-1 system, poses severe problems due 
to the lack of suitable operators in the state preparation. 
Indeed, although any superposition of two states are al- 
lowed in quantum mechanics, it might be difficult to pre- 
pare this superposition in practice. E.g., if we want to 
superpose two different eigenstates of the same Hamil- 
tonian, one has to act on the states with an operator 
which might not be available in the given physical sys- 
tem; this situation is actually encountered in the spin-1 
system. However, emulating a spin-1 system through the 
triplet sector of two qubits provides a viable alternative: 
the two-qubit operation combined with single-qubit ro- 
tations provide sufficient degrees of freedom to carry out 
all required operations for counting, the state preparation 
and the inverse Fourier transformation. Further emula- 
tion of qudits with qubits, however, should not be done 
in the spin-d sector, as this is a waist of resources. In- 
stead, a straightforward sequence of neighboring energy 
states will do. An important element to realize is that 
the Fourier transformation, which we can handle semi- 
classically 'between' the qudits, has to be fully quantum 
'within' the qudits. Hence the larger the chosen count- 
ing basis d is, the larger is the part of the inverse Fourier 
transformation which is done fully quantum. 

We have demonstrated, that the substitution of the 
semi-classical inverse Fourier transformation for the fully- 
quantum does not entail any disadvantage for the count- 
ing algorithm; this is an important result, both with re- 
spect to the stability of the counting algorithm against 
systematic errors (non-integer counting) as well as its 
application in metrology. It thus appears that the con- 
ditional operation in the quantum Fourier transform can 
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be fully substituted by a measurement combined with 
a conditional operation in the semi-classical scheme. In 
both cases, a non-integer reading may affect some of the 
last digits of the seeked number, but the leading digits 
are not compromized. Furthermore, adding additional 
qubits (digits) allows to trade an extension of the hard- 
ware for a better precision in the output. 

In terms of applications of our counter, we have gener- 
alized the scheme producing multi-particle entanglement 
in a Mach-Zehnder interferometer and have proposed its 
use as a quantum voltage-detector, a particular example 
of an analog-digital converter. We may speculate that 
our quantum counting scheme can be generalized to other 
broadband measurement algorithms and thus contribute 
to other applications in quantum metrology^ 



with normalization {4>\4>) = (Xa|Xa) + (Xb\Xb) = 1- The 
state \4>) is not entangled if and only if |x a ) = c\Xb) or 
one of the states \x a )i \Xb) vanishes. Acting with \J V on 



= | t)ie-^ 2)/4 |Xa) 2 



(A4) 



i)ie ^« 2) /4| Xb)2j 



we find that LL 



<Pb - ip a) see Fig. 



disentangles \<f>) if 9 a = 9\> = 9 and 
(here, the angles 6* 0j b and ip a ^ 
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denote the directions of the second spin described with 
the states \Xa,b))- Furthermore, we need the property 
of U v that it transforms a product state |?/>)i|x)2 into a 
product state if either or |x) 2 is directed along the 
z-axis. 



We thank Klaus Ensslin, Alexei Kitaev, Gerald Mil- 
burn, Renato Renner, and Peter Shor for discussions 
and acknowledge financial support by the Pauli-Center 
at ETH Zurich and the Russian Foundation for Basic 
Research under grant No. 08-02-00767-a. 





Appendix A: Derivation of U p and M for the 
qubit-emulated qutrit 



Motivated by our general discussion in section |TV| we 
start out with the counting basis in the form [we choose 
a balanced state \^q) and apply the counting operator 



l*o) 
l*i) 



U) + ut» 



Ci in Eq. (58) to obtain, up to an overal phase] 
1 

1 

7s 
i 

7s 



l*2> 
l*3> 



f(ltt) 

r(ltt) 
f(ltt) 



(Al) 



27ri/3| 



47ri/3| 



U) 
U) 



=47ri/3| 



= 2tti/3| 



II}- 



These counting states are entangled, whereas those defin- 
ing the computational basis, |0) = |tt); |1) — ltl)j an d 
1 2} = |it)) are not- This feature can be conveniently ex- 
ploited in finding the operators U p and M for preparation 
and readout. 

We begin with the preparation step: rather then find- 
ing Up, we search for the inverse operator 1 which 
disentangles the state |^o)- In order to accomplish this 
task, we have the operators a { ^\ a ( z 1} , <jf\ and the 
two-qubit operator a'^a^ at our disposal. In our con- 
struction below, we will make heavy use of the condi- 
tional rotation 



expH^Vf/4). 



(A2) 



This operator will generate the decisive step in the disen- 
tanglement of the qubit states. Let us consider a general 
(entangled) two-qubit state 



t)l|Xa)2 + |l)l|X&>2 



(A3) 



FIG. 12: Action of U^, in disentangling a state. The two com- 
ponents \xa,b) to be aligned point along the same polar angle 
0. The conditional rotation by ip aligns the two components 
if ip matches the azimuthal angle difference, tp = tpt — ip a . 

In order to find the operator U" 1 disentangling l^o); 
we write the latter in the form 



l*o> = I t) 



It)a + U)a , U>i|t>a 



V3 



+ 



V3 



(A5) 



hence | X a) 2 - (| t>2 + I ih)/VS and \ Xb ) 2 = I t) 2 /V3. 
The sequence of operations shown in Fig. |T3| disentangles 
the state |^o) by aligning |xa)2 and |x&) 2 along the di- 
agonal in the xy-plane and produces the state \ip)i\ t} 2 
with 



(A6) 



The remaining one-qubit operations exp(— ina^ /4) and 
exp(— i6a { ^ ) /2) with 9 = 2 arctan(l/v / 2) then produces 
the computational state | ft) up to a phase, U p | tt) = 
exp(i7r/4)|^o)- One easily verifies that the above se- 



quence of operations produces the operator U p , Eq. (60) 



The construction of the inverse quantum Fourier trans- 
formation M for the readout follows the same scheme, 
i.e., we look for the operator M which disentangles the 
counting states and maps |* 2 ) to |2) = 
to |1) = |t|), and |* ) to |0) = |tt); this will be done in 
three consecutive steps, 



U UiU 2 , 



(A7) 



with Ui, i = 0, 1, 2, appropriate unitary operators serving 
to disentangle the three counting states j = 0, 1, 2, 
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FIG. 13: Rotations generating the operator Up 1 ; shown are 
the rotations acting on the components \x<i)2 oc (| t)2 + | D2) 
and \Xb}2 oc | f)2 in !*o), cf. Eq. \A5\ (the arrows reflect the 
polarization angles of the spin states) . The first rotation (a) 
— > (b) by — 7r/2 around the z-axis rotates \t into the direction 
along the j/-axis. The subsequent two-qubit rotation U„/2, cf. 
(b) — > (c), aligns both components along the a;y-diagonal; this 
step disentangles the state. The next rotation by 7r/4 around 
the z-axis makes the two states \ a ,b point along the j/-axis 
and the subsequent rotation by 7r/2 around the :r-axis aligns 
them parallel to the z-axis. The final two rotations of the 
qubit 1 (not shown) transforms the state t)2 into the 

state e~ l7T ^ 4 \ t)i| t)2, which coincides (up to a phase) with 
the computational state |0) = ft}- 



and producing the simple computational states \j), j — 
0,1,2. 

We start with the disentanglement of ^2), which we 
write in the form 



I*-; 



t>2 + e^ /3 |l>2 , e^ /3 |l)i|t) 2 



V3 



V3 



(A8) 



and require U2 to disentangle ^2) while leaving the 
fourth state 1*3) unchanged (up to a phase). This task 
is accomplished by the operator 



U 2 = 



>/4 p -i™™/8 - 



(A9) 
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where the first expression is used in an implementation 
through single and two-qubit gates and the second pro- 
vides the simpler overall unitary matrix, The first two 
rotations prepare the states \xa)2 and |Xf>)2 to allow the 
two-qubit operator 1)^/2 to align them and hence dis- 
entangle the state, see Fig. [Ml the remaining one-qubit 
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FIG. 14: Rotations generating the operator U2; shown are 
the operations acting on the components \Xa)2 oc (| t}2 + 
e 14 " 73 ! D2) and \ X bh « I t>2 in |* 2 ) of qubit 2 in Eq. (til 
(vectors mark the polarizations of the spin states) . The first 
rotation (a) — > (b) by 2-7r/3 around the z-axis rotates the in- 
plane component \a into the x-axis. The second rotation by 
— 7r/2 around the rr-axis makes \b point along the y-axis, cf. 
(b) — > (c). The subsequent two-qubit rotation U 7r /2 , cf- (c) — > 
(d), aligns both components along the x?/-diagonal; this step 
disentangles the state. The remaining rotations of qubit 2 by 
7r/4 around the z-axis and by 7r/2 around the x-axis serve to 
prepare the operation Ui and take the spin 2 back to the z- 
axis. The state ^3} merely picks up a phase, while the states 
l^o) and remain entangled, cf. Eqs. (All I. 



operators acting on qubit 2 serve to prepare the state 
for the action of Ui. The operator U2 leaves the state 
l*3> = I 44-) parallel to itself, U 2 | 44-) = exp(wr/3)| \\), 
and transforms \^f 2 ) into the product state 



^ 2) J^-"^)^f"\MUh, (A10) 



The other two states remain entangled, 



U 2 |*o) = 



U 2 |*i) = ( 



|t)i , e-" /3 U) 



V6 Vs 



It), 



(All) 



^ /2 |t)i|l) 2 
V2 



V6 \/3 



z57r/6 I 



t)iU>s 



V2 



Note that the state of qubit 1 is identical in the above 
two wave functions (pull out a factor e ~ 2 W3 m U 2 |<I' 1 )). 
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while leaving the product states in product states, 



UiU 2 |f 2 ) = e 



UiU 2 |* 3 ) 



U)l|t>2, 

;)iU)a. 



(A14) 



The action of the (conditional) rotations in Ui leading to 
the disentanglement of j^o) and j^i) is shown in Fig. 15 



FIG. 15: Rotations generating the operator Ui; shown are 
the rotations acting on the components |xa)i oc (| t)i/v / 2 + 
e"^ 3 ! 1)0 and \ Xb )i oc | t)i in U 2 |*i) and U 2 |* ) of qubit 
1, cf. Eq. ( |A11 1 (vectors mark the polarizations of the spin 
states). The first rotation (a) — > (b) by — n/6 around the 
z-axis rotates the component \xa) into the yz-plane. The sec- 
ond rotation by —8 around the i-axis takes the two states 
symmetrically around the z-axis, cf. (b) — > (c). The subse- 
quent two-qubit rotation U-n-, cf. (c) — ¥ (d), aligns both com- 
ponents in the xz-plane; this step disentangles the state. The 
remaining rotations of qubit 1 serve to align its state along 
the z-axis. 



Next, the operator Ui, 
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where 9 ~ arctan(V2), takes the remaining two entan- 
gled states U2 1 ^0} and U2 1 ^l) into product states, 
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Finally, the operator Uo, 
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V2 
\ V2 J 
acts as a Hadamard operation on the qubit 2 if the qubit 
1 is in the | f) 1 state and leaves the qubit 2 unchanged if 
the qubit 1 is in the state | hence Uo is a controlled 
Hadamard. The combined action M of the operators U 2 , 
Ui, and Uo finally take the counting states into the com- 
putational basis states, 
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which is nothing but the desired inverse Fourier transfor- 
mation (up to phases). 
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